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Symplectic  and  multisymplectic  methods  for  the  KdV  equation? 

by 

Uri  M.  Ascher 

Department  of  Computer  Science 
University  of  British  Columbia 
Vancouver,  B.C.  V6T  1Z4 
Canada 

ascher@cs.ubc.ea 


Abstract 


The  design  and  development  of  symplectic  methods  for  Hamiltonian  ordinary  differential  equa¬ 
tions  (ODEs)  has  often  yielded  very  powerful  numerical  schemes  with  beautiful  geometric  prop¬ 
erties.  Symplectic  and  other  symmetric  methods  have  been  noted  for  their  superior  performance, 
especially  for  long  time  integration.  On  the  other  hand,  it  has  long  been  known  that  conservative 
discretization  methods  for  nonlinear,  nondissipative  partial  differential  equations  (PDEs)  governing 
wave  phenomena  tend  to  become  numerically  unstable,  especially  for  long  time  integration.  Thus, 
the  common  beliefs  of  two  established  communities  seem  headed  to  a  clash  upon  examining  sym¬ 
plectic  time  discretizations  for  Hamiltonian  semi-discretizations  and  multisymplectic  discretizations 
of  certain  Hamiltonian  PDEs. 

We  therefore  examine  some  symplectic  and  multisymplectic  methods  for  the  notorious  KdV 
equation  with  the  question  whether  the  added  structure  preservation  that  these  methods  offer  (at 
a  price)  is  key  in  providing  higher  quality  methods  for  such  problems.  Please  stay  tuned  for  the 
answer. 
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High  order  A-stable  numerical  methods  for  stiff  problems 

by 

John  Butcher 

Department  of  Mathematics 
The  University  of  Auckland 
Private  Bag  92019,  Auckland 
New  Zealand 

butcher  @math  .auckland .  ac.  nz 


Abstract 


The  only  way  to  obtain  high  order  A-stable  methods  is  by  having  many  stages  of  calculation  in 
each  time  step.  Let  A  denote  the  coefficient  matrix  relating  the  stage  values  to  the  stage  derivatives. 
To  lower  implementation  costs,  specific  structures  for  A  have  to  be  insisted  on.  We  consider  singly- 
implicit  Runge-Kutta  methods  for  which  <r(A)  =  {A},  with  A  chosen  by  stability  criteria  and  general 
linear  methods  which  possess  the  stronger  property  that  A  is  lower  triangular  with  constant  A  on 
the  diagonals.  This  talk  is  an  overview  of  these  two  types  of  method  with  special  emphasis,  in  the 
Runge-Kutta  case,  to  effective  order  and,  in  the  general  linear  case,  to  inherent  RK  stability. 
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Fourth-Order  ESDIRK  schemes  for  Fluid  Mechanics  Applications 

by 

Mark  Carpenter 
NASA  Langley 
carpentr@tabOO.larc.nasa.gov 


Abstract 


Recent  experiences  with  Explicit  Singly  Diagonal  Implicit  Runge-Kutta  (ESDIRK)  schemes  are 
presented  in  the  context  of  high  Reynolds  number,  turbulent  flows.  Structured  and  unstructured 
codes  are  used  to  compare  the  ESDIRK  schemes  with  the  BDF2,  BDF3  and  MEBDF4  schemes. 
Cost  comparisons  are  made  at  different  precisions,  and  the  effects  of  order-reduction  is  quantified. 
Critical  areas  are  identified  for  future  work. 
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Preconditioning  Strategies  for  Models  of  Incompressible  Flow 

by 

Howard  Elman 
University  of  Maryland 
Department  of  Computer  Science 
College  Park,  MD  20742 
elman@cs.umd.edu 


Abstract 


We  describe  some  new  preconditioning  strategies  for  handling  the  algebraic  systems  of  equations 
that  arise  from  discretization  of  the  incompressible  Navier-Stokes  equations.  We  show  how  these 
approaches  relate  to  classical  techniques  such  as  projection  methods  and  SIMPLE,  demonstrate 
how  they  adapt  in  a  straightforward  manner  to  decisions  on  implicit  or  explicit  time  discretization, 
and  explore  their  use  on  a  collection  of  benchmark  problems. 
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SSP  Runge-Kutta  Time  Discretizations 

by 

Sigal  Gottlieb 

University  of  Mass  -  Dartmouth 
&  Brown  University 
Division  of  Applied  Math 
Providence,  RI  02912 
sg@cfm.brown.edu 


Abstract 

In  this  talk  I  will  review  and  and  describe  the  ongoing  study  of  a  class  of  strong  stability  pre¬ 
serving  (SSP)  high  order  time  discretizations  for  semi-discrete  method  of  lines  approximations  of 
partial  differential  equations.  Previously  called  TVD  (total  variation  diminishing)  time  discretiza¬ 
tions,  these  high  order  time  discretization  methods  preserve  the  strong  stability  properties  of  first 
order  Euler  time  stepping  and  has  proved  very  useful  especially  in  solving  hyperbolic  partial  dif¬ 
ferential  equations.  In  this  talk  I  will  describe  the  construction  of  optimal  explicit  SSP  linear 
Runge-Kutta  methods,  the  systematic  study  of  explicit  SSP  Runge-Kutta  methods  for  nonlinear 
problems,  and  the  study  of  the  strong  stability  preserving  property  of  implicit  Runge-Kutta  and 
multi-step  methods,  and  recent  developments  in  the  field  of  SSP  methods. 
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High  order  time  compact  difference  methods  for  wave 

propagation  problems 

by 

Bertil  Gustafsson 
Dept  of  Scientific  Computing 
Uppsala  University 
Uppsala,  Sweden 
bertil@tdb.uu.se 


Abstract 


It  is  well  known  that  wave  propagation  problems  over  long  time  intervals  require  high  order 
methods,  and  many  such  methods  have  been  constructed.  However,  to  keep  high  order  accuracy  also 
in  time  for  problems  with  variable  coefficients  and  minimal  storage,  the  methods  quickly  become 
complicated.  In  this  talk,  we  will  consider  acoustic  wave  propagation,  and  discuss  a  few  finite 
difference  methods  for  the  wave  equation  in  scalar  form  as  well  as  in  first  order  system  form. 

One  way  to  achieve  high  accuracy  in  time,  is  to  use  Taylor  expansion  and  substitute  time 
derivatives  by  space  derivatives  by  using  the  differential  equation.  It  turns  out  that  the  first  order 
system  form  leads  to  less  complicated  approximations,  while  still  keeping  the  conservation  properties 
of  the  continuous  formulation.  We  will  show  that  the  method  works  well  even  for  discontinuous 
coefficients  without  any  special  procedures  across  the  material  interfaces. 

The  deferred  correction  principle  is  well  known  for  obtaining  high  order  accuracy  in  space,  but 
it  can  also  be  used  in  the  time  direction.  We  will  present  some  new  results  for  this  type  of  methods 
based  on  implicit  time  discretization  of  first  order  PDE  systems.  By  slightly  weakening  the  stability 
concept,  it  is  possible  to  keep  the  unconditional  stability  even  for  the  higher  order  methods. 
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Important  aspects  of  geometric  numerical  integration 

by 

Ernst  Hairer 

Section  de  mathematiques 
2-4,  rue  du  Lievre 
CH-1211  Geneve  24 
Switzerland 

Ernst. Hairer@math.unige.ch 


Abstract 


A  short  review  about  properties  and  advantages  of  symplectic,  symmetric  and  reversible  meth¬ 
ods  for  the  numerical  integration  of  Hamiltonian  systems  will  be  given.  The  emphasis  is  on  aspects 
that  are  important  but  sometimes  overlooked,  such  as:  existence  of  globally  defined  modified  Hamil¬ 
tonian,  limits  of  backward  error  analysis,  use  of  variable  step  sizes.  Numerical  phenomena  will  be 
reported  that  still  miss  a  theoretical  explanation. 
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Construction  and  Implementation  of  General  Linear  Methods  for 

Ordinary  Differential  Equations 

by 

Zdzislaw  Jackiewicz 
Department  of  Mathematics 
Arizona  State  University 
Tempe,  Arizona  85287 
jackiewicz@asu.edu 


Abstract 


In  the  first  part  of  this  lecture  we  will  give  the  overview  of  different  approaches  to  the  construc¬ 
tion  of  diagonally  implicit  multistage  integration  methods  for  both  nonstiff  and  stiff  differential 
systems  of  ordinary  differential  equations.  The  identification  of  high  order  methods  with  appro¬ 
priate  stability  properties  requires  the  solution  of  large  systems  of  nonlinear  equations  for  the 
coefficients  of  the  methods.  For  low  orders  these  systems  can  be  generated  and  solved  by  symbolic 
manipulation  packages.  For  high  orders  the  approach  to  the  construction  of  such  methods  is  based 
on  the  computation  of  the  coefficients  of  the  stability  function  by  a  variant  of  the  Fourier  series 
method  and  then  solving  the  resulting  large  systems  of  polynomial  equations  of  high  degree  by  least 
squares  minimization.  Using  these  approaches  both  explicit  and  implicit  methods  were  constructed 
up  to  the  order  eight  with  good  stability  properties  (Runge-Kutta  stability  for  explicit  methods, 
A-stability  and  L-stability  for  implicit  methods). 

In  the  second  part  of  this  talk  we  will  address  different  issues  related  to  the  implementation  of 
general  linear  methods.  They  include  selection  of  initial  stepsize  and  starting  values,  computation 
of  Nordsieck  representation,  efficient  and  reliable  estimation  of  the  local  discretization  errors  for 
nonstiff  and  stiff  equations,  step  size  ond  order  changing  strategies,  construction  of  continuous 
interpolants,  and  updating  vector  of  external  approximations  to  the  solution.  Experiments  with 
variable  step  variable  order  experimental  Matlab  codes  for  both  nonstiff  and  stiff  differential  systems 
on  interesting  test  problems  will  be  presented  and  compared  with  appropriate  codes  from  Matlab 
ODE  suite.  These  experiments  demonstrate  the  high  potential  of  diagonally  implicit  multistage 
integration  methods,  especially  for  stiff  systems  of  differential  equations. 
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Time-Integration  Methods  in  Computational  Aerodynamics 

by 

Antony  Jameson 

Department  of  Aeronautics  and  Astronautics 
Stanford  University 
Stanford,  CA 

j  ameson@baboon  .Stanford .  edu 


Abstract 
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Semi-Lagrangian  Methods  for  high  Reynolds  number  flows 

by 


George  Em  Karniadakis 


Division  of  Applied  Mathematics 
Brown  University 
182  George  Street 
Providence,  RI  02912 
gk@dam.brown.edu 


Abstract 


A  review  of  semi-Lagrangian  methods  for  advection-diffusion  and  incompressible  Navier-Stokes 
equations  will  be  presented.  The  focus  is  on  constructing  stable  schemes  of  second-order  temporal 
accuracy,  as  this  is  a  crucial  element  for  the  successful  application  of  semi-Lagrangian  methods 
to  turbulence  simulations.  Theoretical  results  show  a  non-monotonic  variation  of  the  error,  and 
numerical  results  validate  this  error  estimate  for  the  advection-diffusion  equation.  Two-  and  three- 
dimensional  laminar  and  transitional  flow  simulations  suggest  that  semi-Lagrangian  schemes  are 
more  efficient  than  their  Eulerian  counterparts  for  high-order  discretization  on  non-uniform  grids. 
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Equation-Free  Multiscale  Computation: 

Enabling  Microscopic  Simulators  to  Perform  System  Level  Tasks 

by 

Ioannis  G.  Kevrekidis 
Department  of  Chemical  Engineering 
PACM  and  Mathematics 
Princeton  University 
Princeton,  NJ  08544 
yannis@princeton.edu 


Abstract 


I  will  present  and  discuss  a  framework  for  computer-aided  multiscale  analysis,  which  enables 
models  at  a  ’’fine”  (microscopic/stochastic)  level  of  description  to  perform  modeling  tasks  at  a 
’’coarse”  (macroscopic,  systems)  level.  These  macroscopic  modeling  tasks,  yielding  information 
over  long  time  and  large  space  scales,  are  accomplished  through  appropriately  initialized  calls  to 
the  microscopic  simulator  for  only  short  times  and  small  spatial  domains:  ’’patches”  in  macroscopic 
phase  space-time. 

Traditional  modeling  approaches  start  by  deriving  macroscopic  evolution  equations  (balances 
closed  through  constitutive  relations).  An  arsenal  of  analytical  and  numerical  techniques  for  the 
efficient  solution  of  such  evolution  equations  (usually  Partial  Differential  Equations,  PDEs)  is  then 
brought  to  bear  on  the  problem.  Our  equation-free  (EF)  approach,  introduced  in  PNAS  (2000) 
when  successful,  can  bypass  the  derivation  of  the  macroscopic  evolution  equations  it  when  these 
equations  conceptually  exist  but  are  not  available  in  closed  form. 

We  discuss  how  the  mathematics-assisted  development  of  a  computational  superstructure  may 
enable  alternative  descriptions  of  the  problem  physics  (e.g.  Lattice  Boltzmann  (LB),  Brownian 
Dynamics  (BD),  kinetic  Monte  Carlo  (KMC)  or  Molecular  Dynamics  (MD)  microscopic  simulators, 
executed  over  relatively  short  time  and  space  scales)  to  perform  systems  level  tasks  (integration 
over  relatively  large  time  and  space  scales,  ’’coarse”  bifurcation  analysis,  but  also  optimization  and 
control  tasks)  directly. 

In  effect,  the  procedure  constitutes  a  systems  identification  based,  ’’closure  on  demand”  com¬ 
putational  toolkit,  bridging  microscopic/stochastic  simulation  with  traditional  continuum  scientific 
computation  and  numerical  analysis.  We  illustrate  these  ’’numerical  enabling  technology”  ideas 
through  examples  from  chemical  kinetics  (LB,  KMC),  rheology  (Brownian  Dynamics),  homog¬ 
enization  and  the  computation  of  ’’coarsely  self-similar”  solutions,  and  discuss  various  features, 
limitations  and  potential  extensions  of  the  approach. 
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Jacobian-Free  Newton-Krylov  methods  for  the  accurate  time 
integration  of  stiff  wave  problems. 

by 

Dana  Knoll 

Los  Alamos  National  Laboratory 
M.S.  D413 

Los  Alamos,  NM  87545 
nol@lanl.gov 


Abstract 

We  present  results  from  our  research  on  using  JFNK  methods  for  the  implicitly  balanced  time 
integration  of  stiff  wave  problems.  We  draw  our  examples  from  the  areas  of  geophysical  fluid 
dynamics  and  magnetohydrodynamics.  Stiff  wave  problems  are  a  classical  multiple  time  scale 
problem  where  the  fast  time  scale  (stiff  wave)  physics  is  in  a  “near  balance”  with  other  normal 
modes  in  the  system  resulting  in  a  slow  dynamical  time  scale.  In  this  situation  one  would  like  to 
use  time  steps  on  the  order  of  the  dynamical  scale  for  time  integration.  Historically,  Semi-Implicit 
(SI)  methods  have  been  developed  to  step  over  the  stiff  wave  time  scale  in  a  stable  fashion.  However, 
SI  methods  require  some  linearization  and  time  splitting,  and  both  of  these  can  produce  additional 
time  integration  errors.  Our  algorithmic  approach  employs  SI  methods  as  preconditioners  to  JFNK, 
resulting  in  an  implicitly  balanced  method  (no  linearization  or  time  splitting).  In  this  presentation 
we  will  provide  an  overview  of  SI  methods  in  a  variety  of  applications.  We  will  present  details  of 
our  algorithmic  approach.  Then  we  will  compare  accuracy  and  efficiency  of  our  method  with  more 
traditional  SI  methods. 
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Time-integration  of  low  Mach  number  reacting  flow  with  detailed 

chemical  kinetics 

.  by 

Habib  Najm 

Sandia  National  Laboritories 
7011  East  Avenue,  MS  9051 
Livermore,  CA  94550 
hnnajm@ca.sandia.gov 


Abstract 

The  time-integration  of  low  Mach  number  reacting  flow  with  detailed  kinetics  presents  a  number 
of  challenges.  Even  laboratory-scale  realizations  of  these  flows  with  simple  fuels  are  characterized 
by  large  ranges  of  length  and  time  scales  as  well  as  significant  complexity  in  chemical  source  term 
evaluations.  Efficient  time  integration  of  these  systems  requires  attention  to  be  paid  to  both  spa¬ 
tial  and  temporal  discretization  constructions.  In  this  talk  I  will  present  our  experiences  in  this 
field,  outlining  primarily  time  integration  issues,  involving  projection-scheme  integration  of  mul¬ 
tidimensional  low  Mach  number  flames,  the  coupling  of  momentum  and  species  time  integration 
procedures,  as  well  as  the  utility  of  semi-implicit  (IMEX)  and  operator-split  time-integration  con¬ 
structions.  Given  the  significant  costs  associated  with  the  time  integration  of  chemical  source  terms, 
I  will  also  discuss  this  area  specifically,  and  will  address  various  means  of  model  reduction  and/or 
solution-tabulation  to  reduce  overall  time  integration  costs. 
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Abstract 


We  present  implicit-explicit  (IMEX)  Runge  Kutta  methods  for  hyperbolic  systems  of  conser¬ 
vation  laws  with  stiff  relaxation  terms.  Here  we  restrict  our  analysis  to  diagonally  implicit  Runge 
Kutta  (DIRK)  schemes  and  consider  schemes  up  to  order  3  that  are  asymptotic  preserving  (AP) 
and  strong-stability-preserving  (SSP)  in  the  zero  relaxation  limit.  Accuracy  and  stability  properties 
of  these  schemes  are  studied  both  analytically  and  numerically.  The  schemes  are  then  combined 
with  a  weighted  essentially  non  oscillatory  (WENO)  finite  difference  strategy  in  space  to  obtain 
fully  discrete  high  order  schemes.  Several  applications  to  hyperbolic  relaxation  systems  in  stiff  and 
non  stiff  regimes  show  the  efficiency  and  the  robustness  of  the  present  approach. 
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Abstract 


This  talk  is  about  the  numerical  solution  of  initial  value  problems  for  systems  of  ordinary 
differential  equations  (ODEs).  At  first  these  problems  were  solved  with  a  fixed  method  and  constant 
step  size,  but  nowadays  the  general-purpose  codes  vary  the  step  size,  and  possibly  the  method,  as 
the  integration  proceeds.  Estimating  and  controlling  some  measure  of  error  by  variation  of  step 
size/method  inspires  some  confidence  in  the  numerical  solution  and  it  makes  possible  the  solution 
of  hard  problems.  Common  ways  of  doing  this  are  explained  briefly  in  the  talk. 
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Abstract 


The  electromagnetic  force  becomes  a  new  dimension  for  enhancing  the  aerodynamic  performance 
of  aerospace  vehicles.  In  order  to  simulate  the  interdisciplinary  phenomena,  the  Navier-Stokes  and 
Maxwell  equations  in  the  time  domain  must  be  integrated  on  the  same  frame  of  reference.  For  a  wide 
range  of  application  from  hypersonic  flow  control  to  microwave  propagation,  the  resultant  nonlinear 
partial  differential  equations  offer  a  formidable  challenge  for  numerical  analysis.  The  experience 
and  physical  insight  using  the  approximate  Riemann  and  Compact-Differencing  formulation  as  well 
as  several  temporal  discritizations  will  be  shared.  The  most  recent  development  and  advancement 
in  numerical  procedures  for  solving  this  system  of  governing  equation  are  delineated. 
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Abstract 


We  describe  the  implementation  of  a  new  dynamical  core  for  atmospheric  general  circulation 
models.  To  avoid  singularities  at  the  poles,  a  cubed-sphere  coordinate  system  is  employed.  A  high- 
order  spectral  element  discretisation  in  the  horizontal  direction  is  combined  with  a  second  order 
finite  difference  scheme  in  the  vertical.  Large  time  steps  beyond  the  advective  CFL  limit  are  made 
possible  by  the  operator-integration  factor  (OIF)  scheme  of  Maday  et  al  (1990).  A  three-time  level 
leapfrog  scheme  for  advection  is  combined  with  a  Crank-Nicholson  (semi-implicit)  scheme  applied 
to  the  terms  responsible  for  gravity  waves.  Linear  stability  analysis  reveals  that  the  resulting 
discretized  system  of  PDE’s  is  neutrally  stable  for  the  meteorologically  important  Rossby  waves. 
Variable  resolution  has  been  implemented  recently  using  non-conforming  spectral  elements  with 
interpolation  at  element  boundaries.  Numerical  results  are  reported  for  the  baroclinic  instability 
test  problem  proposed  by  Polvani  et  al  (2003). 
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Abstract 


For  solving  highly  stiff  diffusion-reaction  systems  one  often  adopts  the  implicit  method  of  lines 
approach.  The  main  issue  then  is  how  to  solve  the  arising  systems  of  linear  and  nonlinear  algebraic 
equations.  In  higher  space  dimension  this  can  be  costly  and  problematical.  In  this  talk  we  avoid  this 
and  adopt  the  explicit  method  of  lines  approach  and  focus  on  stabilized  methods  from  the  Runge- 
Kutta-Chebyshev  family.  An  implicit-explicit  (IMEX)  extension  of  an  existing  explicit  Runge- 
Kutta-Chebyshev  (RKC)  scheme  designed  for  parabolic  PDEs  is  proposed  for  diffusion-reaction 
problems  with  severely  stiff  reaction  terms.  The  IMEX  scheme  treats  these  reaction  terms  implicitly 
and  diffusion  terms  explicitly.  Within  the  setting  of  linear  stability  theory,  the  new  IMEX  scheme 
is  unconditionally  stable  for  reaction  terms  having  a  Jacobian  matrix  with  a  real  spectrum.  For 
diffusion  terms  the  stability  characteristics  remain  unchanged.  A  numerical  comparison  for  a  highly 
stiff,  nonlinear  radiation-diffusion  problem  between  an  RKC  solver,  an  IMEX-RKC  solver,  and  the 
implicit  BDF  solver  VODPK  using  the  Krylov  solver  GMRES  illustrates  the  excellent  performance 
of  the  new  IMEX  scheme. 
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